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Numerical study of the scaling of transmission fluctuations in the 1-D localization problem in 
the presence of absorption is carried out. Violations of single parameter scaling for lossy systems are 
found and explained on the basis of a new criterion for different types of scaling behavior derived 
^ 1 by Deych et al [Phys. Rev. Lett., 84, 2678 (2000)]. 
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I. INTRODUCTION 



The single parameter scaling (SPS) hypothesis is the cornerstone of the current understanding of the localizatio 
phenomena. It was originally formulated in terms of the scaling behavior of the conductance of disordered conductors,! 
where it was suggested that when the length, L, of a disordered conductor increases, the evolution of the conductance, 
g(L), is determined by a single parameter, g itself. For one-dimensional systems Landauer's formulao expresses 
C/3 \ the electron conductance in terms of electron reflection, R{L) = rr* , and transmission, T(L) = tt*, coefficients 
as g(L) = R(L)/T(L), and thereby allows the considering of electron transport on an equal footing with these 
propagation of, for example, light. It was recognized later that SPS must be understood in terms of properties of the 
C$ ■ entire distribution of g (or T), and that the most appropriate quantity to deal with isa 
P ' 

I i m = (1/2L) ln(l + l/g{L)) = (1/2L) ln(l/T(L)). 

■ 

In the limit of large L, this parameter is normally distributed with the average 7 = (7(£)) = \miL^ 00 ; y(L) and.-the 
variance c 2 (L) = (7 2 (L)) — j 2 . The limiting value of j(L) is known in the theory of products of random matricesQ as 
i^i - the Lyapunov exponent (LE). The inverse quantity, the localization length, l[ oc = 7 _1 , determines the main length 
scale in the localization regime.! SPS in this context means that the variance of 7 is not an independent parameter, 
■ but it is determined by 7 itself implying a simple relationship between the twcB 
> ■ 

<r 2 = 7 /i. (1) 

This expression was obtained for the one-dimensional model in Ref. |^ assuming complete randomization of the phases 
of complex transmission and reflection coefficients over a microscopic length scale l p h <C hoc (phase randomization 
I 1 hypothesis). Later the phase randomization hypothesis was used by many different authors to rederive Eq. ([!]) (see, 
for instance, Refs. |5|||), and the inequality l p h <§C h oc came to be regarded as the criterion for SPS. However, there 
"jS ■ were earlier signs that the phase randomization hypothesis is neither a necessary nor a sufficient condition for SPS to 
occur. For instance, numerical simulations of Ref. |6| and analytical calculations of Ref. [?] showed that in the center of 
a conductivity band of the 1-1? Anderson model, SPS holds even though the phase is not randomized, provided that 
the disorder is weak. Later, numerical simulations of a random, periodic-on-average modelu demonstrated a strong 
violation of SPS in the band-gaps of the spectrum of the underlined system without disorder, which existed even for 
weak disorder and, actually, diminished, when disorder increased. This was contrary to the behavior found for states 
from the original conducting band, for which minor deviations of the variance from Eq. ([l]) occurcd when disorder 
r* ■ becomes strong enough.Q 

The final realization of the fact that the phase randomization hypothesis has nothing to do with SPS came in Ref. 
|9|. In that paper, the variance of LE was calculated exactly for the Lloyd model and the SPS equation (|l|) was derived 
without ad hoc assumptions. It was found that the emergence of SPS is governed by a new length scale / s , related to 
the integral density of states with the new criterion for SPS being k — h oc /h ^S> 1< On the basis of the exact solutions 
it was conjecturedu that in the region of the spectrum close to its original boundary, the parameter l s can be defined 
in a generic case as 

l s = a/N(E), (2) 

where a is the lattice constant and N(E) is the number of states between the lowest genuine boundary of the spectrum 
of the disordered system and E, normalized by the total number of states in the band (such that < N(E) < 1). 
Numerical studies undertaken in Ref. |^ evidenced that this assumption was valid for the Anderson model and the 
periodic-on-average model with rectangular distribution of random parameters. An additional implication of the 
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results of Ref . || is that the random matrix theory approach, which also reproduces Eq. ([j]) ,E30 does not apply to 
spectral regions with depleted differential density of states. 

The objective of this paper is to show that two different scaling regimes governed by the parameter l s found in Ref. ^ 
exist also in disordered systems with absorption. Having in mind applications to light propagation in random photonic 
band gap materials, we consider how the inclusion of small absorption affects scaling properties of transmission. We 
show that the conjectured definition of l s in terms of the integral density of states can also be applied to disordered 
systems with absorption. 

Within the phase randomization hypothesis approach, statistics of the transmission in lossy one-dimensional di- 
electrics was considered analytically in a number of papers.tjO Following Ref. |l2] the relation between the variance 
and the localization length in the presence of absorption can be presented in the form 

r = r (/3) = l + 2/3e 2/3 K(-2/3), (3) 

where r = a 2 L/j, (3 = li oc /la, and l a is the absorption length in the absence of randomness, 

Ei{x) = [ dt exp(t)/t 



is the exponential integral. In the absence of absorption, Eq. (||) reduces to the regular SPS form r = 1. In two 
limiting cases Eq. (||) gives the asymptotes r = 1 — 2/31n(l//3), (3 <C 1 and r = 1/2/3, /3 >■ 1. 

Presuming that the length scale l s retains its meaning in the case under consideration, we expect that the parameter 
r deviates from the phase randomization hypothesis prediction to(/3) [Eq. @] in the vicinity of the boundaries of 
the spectrum in accordance with the same criterion li oc 3> l s as in Ref. 0. Using numerical simulations of a periodic- 
on-average one-dimensional lossy system, we show that, indeed, the parameter k = li oc /l s sets a valid criterion for 
validity of Eq. (^). As an additional benefit, we demonstrate that within the range of its validity, Eq. represents 
a universal, model independent relation between the variance and the mean value of LE. Though the statistics of 
transmission is determined now by two parameters, one can still regard Eq. (||) as a generalization of SPS because 
the absorption length is independent of disorder. The deviations from the phase-randomization based results of Ref. 
|l2| studied in our paper must be clearly distinguished from results of Ref. [14]. In the latter paper brake down of the 
phase randomization was obtained in the case of very strong disorder and strong absorption for states at the center 
of the original band. The results of our paper indicate that (i) violation of the generalized single parameter scaling 
occurs at weak disorders for states close to the band edge of the original spectrum, (ii) this violation is not related to 
the phase randomization but is controlled by the parameter l s . 



II. MODEL AND THE METHOD OF CALCULATIONS 



We consider a classical transverse electromagnetic wave propagating normally through a stack of alternating dielec- 
tric slabs with dielectric constants ei and £2- The widths of the stacks of the first kind is distributed uniformly in the 
interval (d\ — 8,di + 5) while the width of the others is being kept constant cfo. The propagation of the waves in the 
superlattice consists of free propagation in the slabs and scattering at the interfaces, where the boundary conditions 
should be satisfied. It can be described using the transfer matrix formalism for the vector v n = (E n , E n /k), where 
E n ,E n are the electric field and its derivative at n-th interface and k = lu/c. The presence of absorption can be 

accounted for by adding a constant complex part to the dielectric functions: e\ — e[ °\l + 
where a is a damping coefficient. Vectors on neighboring interfaces are connected via the transfer matrix 



and €2 = e 2 (1 + ict), 



T _ ( cos k n d n 
\ &n sin k n d n 



jr" sin k n d n 
cos kj^dji 



(4) 



where k n = k^JT^. The transfer matrix of the entire system is T(a, L) — II^ 1 T„(a), where L = N(d\ + cfe). LE is 
defined through the transmission coefficient for the superlattice: 



7(L,a) = -— (ln(tt*)) = — ( In 



2detT(a,L) 



(Tu(a, L) + T 22 (a, L)) - i(Ti 2 (a, L) - T 21 (ar, L)) 



(5) 



here (...) denotes the average over an ensemble of configurations. We find that using this definition in numerical 
simulations has one significant shortcoming. In long systems , the transmission coefficient tt* falls bellow the computer 
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roundoff accuracy. The usual remedy for this problem in the absence of absorption is to use an alternative definition 
of 7 



1 / \\T(L)v \\ 2 
7(L ' aH 4AF\ ln T^ 



(6) 



which allows one to consider very long systemsfl In the presence of absorption, this definition must be generalized 
because the simple substitution of the transfer-matrix, Eq. (||), in this equation would lead to a wrong result. The 
problem is that an exponentially growing eigenvalue of the transfer matrix appears in the denominator of Eq. (|^), 
while in Eq. ([jj) it is in the numerator. Therefore, the contribution from absorption enters the final answer for the 
LE with different relative signs in these two equation. We argue that Eq. (^) must be modified in order to agree with 
the original definition of 7 as follows: 



j(L,a) 
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where vq is a vector of general position. The proof of equivalency of Eq. 
diagonalize the transfer matrix T(a,_L): 



(7) 



and Eq. (M) goes as follows. First, we 



T(a, L) = T {D \a,L) U = 



e v 2 {a)2N 



U, 



(8) 



here U and are some unitary matrices. For argument sake, we assume that |^i(a)| > |f2(a)|. Next, we notice the 
following relation between eigenvalues of the transfer matrix 



vi{a) = -v<2 (-a), 



(9) 



which follows simply from reciprocity of our system. Indeed, propagation of the waves in opposite direction should 
be described by the matrix T (a,L) = [T(— a, Z)] -1 , that leads to Eq (||). Now, substituting (a, L) into Eqs. 
(||) and (0) (corrections due to matrices U, W are negligible in the limit TV — > 00) one obtains j(L, a) = —V2(a) and 
"f(L,a) = isi(a) respectively. Taking into account the reciprocity relation Eq. (|9|) we establishes the equivalency of 
Eqs. (§) and P). 



III. RESULTS 




P 

FIG. 1. Open circles depict the parameter r(/3) computed for the periodic-on-average system of 100,000 layers far from the 
band edge. 0.86 < k < 1.4, the disorder parameter is <5 = 0.45, the damping a = 0.00125. The solid line shows To(/3) given by 
Eq. i) 
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In numerical simulations we used the following set of parameters e\ = 1.2, 62 = 1, (di) = = 1. The disorder 
parameter, (5, and the absorption rate, a, were variable parameters. To calculate the moments of 7, we averaged the 
characteristics of systems as long as 100, 000 layers over 5000 realizations. The size of the stack was chosen to be at 
least five times the localization length or the absorption length. In the ordered system the first forbidden gap lies 
between k — 1.456 and k — 1.543. Since the localization length depends on the frequency of the wave it is possible to 
study the function t(/3(oj)) by changing the frequency. First, we compare our numerical results with the analytical 
formula of Ref. |l2| in the region of frequencies well inside the first allowed band, where we expect these results to 
coincide. Fig. 1. shows excellent agreement between the computed t(/3) for 0.86 < k < 1.4 and Eq. (||). 




FIG. 2. The dependence of the lyapunov exponent (dashed lines) and its variance (solid lines) averaged over 2500 realizations 
with and without absorption as function of the frequency on the band edge. The length of the system 10,000, disorder parameter 
5 = 0.25, damping 7 = 0.0025. 

The localization length decreases rapidly when the frequency approaches the band edge, while LE, consequently, 
increases. As follows from Eq. (|l|), in the absence of absorption, a 1 should follow the LE. Fig. 2 depicts the dependence 
of LE and its variance with and without absorption. As one can see, at the band edge, k — 1.456, SPS breaks down. 
Indeed, while LE grows with the increasing frequency, the variance drops in both cases, with and without absorption. 




FIG. 3. t(0) (crosses), r(/3) (open circles) and to(/3) (solid line) on the band edge plotted as function of frequency for the 
same set of parameters as on Fig. 2 
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In order to compare numerical and analytical results in the presence of absorption, we calculated the parameter r as 
a function of frequency k, using the data presented in Fig. 2. The results are presented in Fig. 3, from which one can 
see that the computed r(/3) deviates from r (/3) that represents Eq. (||), and that this deviation occurs at the same 
frequency at which r(0) deviates from unity. These graphs convincingly demonstrate, that even in the presence of 
absorption, the spectrum of the system is separated in groups with different scaling properties, and that the boundary 
between the groups coincide with the boundary of the original spectrum. 

The next question we need to address is whether the transition between the scaling regimes is governed by the same 
parameter l s that was introduced in Ref. ^. According to that paper, the new parameter l s becomes greater than the 
localization length at the band edge, which results in a deviation from SPS. Using definition of l s suggested in Ref. || 
regarding the integrated density of states, Eq. (||), we numerically calculated this parameter for the system studied 
in the present paper. The density of states was calculated using the phase formalism (see, for instance, Ref. for 
the system without absorption. Fig. 4 shows all relevant length parameters: localization length, absorption length 
and l s . 
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k 

FIG. 4. Length scales for the system with parameters used on Fig. 2,3. Solid line depicts l a , dashed line is the localization 
length hoc, and dotted line is the absorption length l a 

At the band edge, l s grows rapidly because very few new states appear within the former bandgap, and N(E) must 
already be close to unity at the band edge. It reaches unity at a new fluctuation boundary of the spectrum near 
the center of the gap. We assume that the disorder is not very strong such that the fluctuation boundaries inside 
former bandgaps do exist. If the disorder is strong enough, or if its statistical properties are such that the entire 
bandgap is filled with fluctuation states a definition of l s is still possible, but the situation becomes more complicated, 
and we do not consider it here. At certain points l s grows larger than the localization length li oc , and one can find 
comparing Figs. (3) and (4), that r(/3) starts deviating from To(/3) at the same frequency. In order to make this more 
transparent, it is convienient to plot various r's versus k — k oc /l s - 

Without absorption, we see the crossover to SPS in its pure form at k ~ 1 reported in Refs. [||]. When absorption is 
present, the crossover still occurs at k ~ 1, but now to the modified SPS behavior to(/3). We stress that the crossover 
occurs while a <C hoc, a <C l a , hoc L, and l a <C L, so that the system remains in the meaningful scaling regime. 
This demonstrates that the condition 1 <C n establishes the criterion for modified SPS even in the presence of the 
absorption. 

To check that our results are not model specific we also studied the Anderson model, with absorption introduced as 
a nonrandom imaginary part of the on-site random energy. In the Anderson model without disorder and absorption, 
there is a single allowed band form E = — 2 to £' = 2. 
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FIG. 5. r(0) (crosses), r(/3) (open circles) and to(/3) (solid line) on the band edge plotted as function of k — li oc /ls for the 
same set of parameters as on Fig. 2,3 

Figure 6 shows that t(/3) deviates from tq(/3) at the edges of the allowed band as would be expected on the basis 
of the results presented for the periodic-on-average system. Inside the conductivity band, numerical and analytical 
results show excellent agreement. This is interesting in itself, since Eq. (|j|) was derived for a continuous model,t2l and 
we see that in the SPS regime it holds also for the tight-binding discrete model. 
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FIG. 6. r(0) (crosses), t(/3) (open circles) and tq((3) (solid line) plotted as function of frequency averaged over 1000 realiza- 
tions. The length of the system 10,000 cells, disorder parameter 5 — 0.25, damping 7 = 0.0025. 

In conclusion, we computed the Lyapunov exponent and its variance for a periodic-on-average layered system and 
for the one-dimensional Anderson model. We studied the deviation from the absorption-modified SPS expression for 
the variance.EJ We showed that the new length scale, introduced in Ref. || in order to explain violation of SPS in 
the systems without absorption, retains its significance when absorption is present. We also showed that the same 
criterion k — li oc jl s 3> 1 derived for lossless systems in Ref. ^ controls scaling behavior in lossy systems. 
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